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Abstract 

The critical behaviour of the one- dimensional q-state Potts model with long-range 
interactions decaying with distance r as r~( 1+<T ) has been studied in the wide range 
of parameters < a < 1 and ^ < q < 64. A transfer matrix has been constructed 
for a truncated range of interactions for integer and continuous q, and finite range 
scaling has been applied. Results for the phase diagram and the correlation length 
critical exponent are presented. 
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1 Introduction 



Recently, finite range scaling (FRS) appeared as an efficient method to deal with 
various discrete one-dimensional models with long range (LR) interactions (Glumac 
and Uzelac, 1989, 1991). The aim of the present study is to apply it to the one- 
dimensional q-state Potts model with LR interactions, defined by the hamiltonian 

where the interaction constant is chosen to be equal to ks and the inverse temper- 
ature is denoted by K. The symbol s$ denotes a Potts variable at site i which can 
assume q different states, while S(si, Sj) is the Kronecker symbol. 

Such a model is expected to have, due to the LR of interactions, a nontrivial 
phase transition at finite temperature, like an Ising model (as its special case q=2), 
although the problem becomes much more complex due to possibility of changing 
degrees of freedom through parameter q. 

The Potts model has been extensively studied in the case of short range (SR) 
interactions. In the two dimensional SR case for example, besides integer-q models 
which are equivalent to the classical (q-l)/2-spin models, there are also non-integer 
q models which are connected with other critical problems such as percolation (q=l 
limit), resistor network (q=0 limit), dilute spin glass (q=l/2 limit), branched poly- 
mers (0<g < 1) (see review by Wu, 1982). 

It would be interesting to extend those concepts to LR interactions. However, in 
contrast to the rich literature on the SR Potts model, the LR Potts model has been 
much less explored. The best explored member of the q-state Potts models is the q=2 
model where both, analytical (Dyson 1969; Fisher, Ma and Nickel, 1972; Kosterlitz, 
1976) and numerical (Nagle and Bonner, 1970; Glumac and Uzelac 1989,1991) results 
are available. Dyson has shown that, as a consequence of the LR forces, there is a 
critical temperature different from zero as long as < a < 1. Fisher et al. have 
performed e = 2a — d and 1/n expansions in the context of the renormalization group 
(RG) for all dimensions d and a ^ 2 on a more general system with an n-component 
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order parameter. They obtained the correlation function exponent r/ = 2 — a to all 
orders in e, v — 1/a when e < (classical regime), v = 1/a + ^§^r + • • ■ when 
e > 0, and leading irrelevant exponent y 3 = — |e|. A similar expansion, in l-a > 0, 
has been done by Kosterlitz (1976), who obtained : 1/V = — y 3 = \J2(1 — a) when 
a — > 1. Very recently, a new method, called cycle expansion, has been demonstrated 
by Mainieri (1992) on the example of the s=l/2 Ising model. Mainieri obtained very 
accurate estimates of the critical temperature and other thermodynamic quantities 
in the critical region. 

The literature on q ^ 2 LR Potts models is relatively poor. By using the RG 
e-expansion, Priest and Lubensky (1976) have found a fixed point and critical ex- 
ponent v to the first order in e = 3a — d. More recently, Theumann and Gusmao 
(1985), performed similar calculation to the second order of e = 3(7 — d. They 
obtained the following results for q< 3: i] = 2 — a to all orders of e, and 



where a = q 2 (q -2),/3 = q 2 (q - 3), 7 = q\q 2 - 6q + 10), S(a) = ^(|<r) - i^(a) - 
ip(a/2) + V'(l)) an d C(°") = r"(cr) • ^ s ^ ne logarithmic derivative of the 

r function. In the q — > 1 limit, their results are applicable to the description of 
percolation critical behaviour for a random Ising ferromagnet with a LR power-law 
interaction. 

The advantage of the FRS method is the capability to cover, unlike the e- 
expansion, a wide range of parameter space. In the present article the phase diagram 
and the correlation length critical exponent will be studied in the region < a < 1 
and jq < q < 64. 

The main idea of our approach is the following: we truncate the originally infinite- 
range interaction to the L first neighbours and solve exactly (although numerically) 
the finite-range version of (1) 

N L 1 

- (3H = £ ^ Ktfisi, s i+j ) - -U c . (3) 
i=i j=i 1 
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where Kj denotes jt+t- Then, by using the FRS method (section 2.1) with the ap- 
propriate extrapolating procedure (section 3), the oo behaviour will be deduced. 

For this purpose a transfer matrix for a reduced model defined by eq. (3) will 
be constructed. Two different formalisms will be applied: one for integer q values 
(section 2.2) and the other for non-integer q values ( section 2.3). 

The results are presented in section 3, while section 4 contains a short discussion 
and some open questions. 

2 Method 
2.1 FRS 

The FRS method has been constructed in analogy with the FSS (Fisher and Barber, 
1972), where instead of finite-size, the finite-range of interactions is scaled (Glumac 
and Uzelac, 1989,1991). The basic idea of FRS is that by studying the sequence of 
systems with their long-range interactions truncated above certain range, one can 
obtain, by using scaling properties, the information on the critical behaviour of the 
true infinite system. 

Let Ax>(t) be the physical quantity of the true long-range system, which alge- 
braically diverges, in the vicinity of the critical point t=0 

Ao.it) ~ Aot- p (4) 

where t=(T-T c )/T c , T c is the critical temperature, p is the related critical exponent, 
and A is a constant. Then, analogous to the FSS hypothesis one can assume that 
for large finite range M and small t, A^W can be written in the form 

A M (t)=A 00 (t).fA (5) 

Soo 

where f is a homogeneous function with following properties 

lim fix) = 1, lim fix) = const. -x^ (6) 

x^oo J w ' x->0 V ' V ' 

By applying the equation (5) to the correlation length £oo(£) = (,ot~ u , the stan- 
dard scaling procedure gives the condition for critical temperature through the fixed 
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M M> 



and the expression for correlation length critical exponent v 




- 1 
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where M'=M-1 in all calculations, and £' is temperature derivative of £. 

There are two important facts concerning the above method which are of interest. 

First, since the critical behaviour is essentially dependent on the range of inter- 
action, one can (unlike in the FSS (Brezin, 1982)) expect applicability of FRS in 
mean-field as well as in the non-trivial region (Glumac and Uzelac 1989, 1991). Sec- 
ond, the correlation length calculated from the transfer matrix presents the average 
distance between domain walls for an infinitely long strip so that relation (7) gives 
the transition temperature for both first and second-order phase transitions. This 
fact is common with FSS (Binder, 1987). 

2.2 Transfer matrix: integer-q formalism 

It is straightforward to construct the transfer matrix for model (3) with an integer 
number of Potts states. 

The chain with a range of interaction L, can be represented as a strip with 
columns of height L (Uzelac and Glumac, 1988). Each column then can be consid- 
ered as an object with q L possible states, interacting only with its first neighbour. 

The transfer matrix T is given by: 



where j m = 0, 1, . . . , q — 1 are the elements of the L-component vector of states \j) 
of a column of height L. The matrix T can be further decomposed into a product 
of L matrices T n , each one describing the addition of one more site to the column 
(Temperley and Lieb, 1971). 





T = T x ..... T L 



(10) 
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L (n-1 L-1 \ 

(i\T n \j) = Yl S(U,ji)exp I K k 5(j n _ k ,j n ) + ^ K k 5(j L+n _ k ,j n ) + K L 5{i ni j n ) \ 

1=1 U=l k=n ) 

(11) 

There is a simple relation between neighbouring one-site matrices: 

U T T„ +1 U = T n (12) 

where T i+1 =T! and U is a matrix of the translation operator in the vertical strip 
direction, with matrix elements given by 

(i\U\j) = 5(i 1 ,j L )5(i 2 ,ji)5(i 3 ,j 2 ), ■ ■ ■ , S(i L - U j L -2)S(i L , 3l-i)- (13) 

The matrix U satisfies the relations U L = 1 and U L_1 = U T = U _1 . Conse- 
quently, the transfer matrix can be written as the L-th power of a single matrix as 
a peculiarity of the present long-range model: 

T = (U • T L ) L = f L (14) 
(i\T\j) = 5(j 1 ,i 2 )5(j 2 ,i 3 ) ■ ..6(J L -i,i L )exp I ^ K L+i-mS(j L ,i m ) \ (15) 

lm=l J 

Notice that ~T has only q nonzero elements per row which greatly reduces computer 
memory. 

Further reduction is obtained by using the symmetry properties of Potts interac- 
tion. In the present calculation, only the translation invariance in the space of Potts 
states has been used, which decomposes ~T into q submatrices (of order q L ~ l ). The 
q-1 of them have the identical eigenvalues as a consequence of Kronecker-<5 type of 
interaction. Thus we have to diagonalize only two submatrices or, more precisely, 
we have to find only the largest eigenvalue of each of the two submatrices, and 
/i 2 ,L respectively, to be able to construct the correlation length £ L 

L 1 



(16) 



where Al > A2 and /il > /i2 are the eigenvalues of T and ~T respectively. 

For this purpose, the method of direct iteration with a single vector is suitable 
(Wilkinson, 1965). If one defines a set of vectors bj in a following way 

Ab, = b, +1 (17) 



for arbitrary initial value bo, where A is a real nonsymmetric matrix, then the 
largest eigenvalue \x of matrix A is given by 

VL = hm -M^) (is) 
j^oc bj-i{n) 

for any vector's component n. By applying this method for cases with q= 2, 3, 4 
and 5, we obtained the two largest eigenvalues of ~T for maximum ranges L= 20, 
13, 10 and 9, respectively. A few ten to a few hundred iterations were required to 
obtain the leading eigenvalue with an accuracy of 10~ 15 for both submatrices. 

2.3 Transfer matrix: continuous-q formalism 

As mentioned in Introduction, there are some interesting cases (q=0, 1/2, 1, ...) to 
which the described integer-q formalism cannot be applied or becomes inefficient 
(limitation to small ranges for large q models). Thus, it would be advantageous 
to formulate a q-independent transfer matrix, which would allow to treat 
continuous parameter. Such a type of transfer matrix has already been constructed 
for the 2D Potts model with SR interaction (Blote, Nightingale, Derrida 1981; Blote 
and Nightingale, 1982). We will shortly explain the application of that method to 
the considered LR model. 

Let us translate the partition function of model (3), into the graph-theory lan- 
guage. First, we rewrite the basic expression exp{Kj • 5(si, s i+ j)} as 

e K 3 8{ Sl ,s l+] ) = l + Vj S( SijSi+j ) : (19) 

where Vj = e Kj — 1, so that the partition function, with periodic boundary conditions, 
becomes 

vee-ee (2°) 

Sl=ls 2 = l S N = lj = l 

[1 + Vj6(si, s 1+j )][l + Vj5(s 2 , s 2+j )] . . . [1 + VjS(s N , s N+j )] 

If we perform the multiplications of the above terms, the one-to-one correspon- 
dence between every term of that expansion and the particular graph on the mxL 
lattice (N=m-L) can be established in the standard way (Blote and Nightingale, 
1982). 



Now, the partition function becomes 

^ = ^-E% l(G) -4 2(G) -..--4 L(G) V (G) . (2i) 

G 

The above summation goes over all possible graphs on the lattice with L different 
types of bonds described by Uj = Vj/q , bj(G) is the number of bonds of type j, and 
1(G) is the number of independent loops on the graph G. 

A certain number of different bonds in every single graph, produce connections 
between sites. If one considers a strip of height L, then all possible interactions 
between all sites in that strip produce a particular state of connectivity a of the sites 
in the last column. The LR interactions produce a larger number of connectivities 
than in the 2D SR models. 

The connectivity is represented by an integer, and it will be used as an index 
of the transfer matrix. The partition function can be written as a vector with 
components Z^(a) , where each component Z^(a) describes the contribution to the 
partition function of all graphs that produce one particular state of connectivity 
a in the last column. By adding a new column, we get Z N+L which again can be 
written as a vector Z N+L ((3). The connection between those two partition functions 
is given by the transfer matrix in the following way 

Z N+L ((3) = J2T((3,a)Z N (a) (22) 

a 

where the summation goes over all graphs, produced by the added column, that 
lead from connectivity a to (5. A transfer matrix decomposition analogous to that 
of section 2.1, leads to 

T = gTi • . . . • qT L = ( g TiV T ) L = T L (23) 

where Tj describes the addition of the j-th site to the column, and V is the matrix 
of the translation operator in the vertical strip direction. Further decomposition of 
Ti is possible 

Ti = T 1>L (u L ) • J] T u(«i) ( 24 ) 

j=L-l 
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Each of the Tij(uj) matrices describes the presence of an Uj type bond on the first 
site and has at most two nonzero elements per row. T 1>l (vll) is an upper triangular 
matrix, while the other matrices are lower triangular. It is easy to construct a matrix 
M which relates the lower triangular matrices 

M T T 1 ,(x)M = T 1 , +1 (x) (25) 

The matrix ~T then becomes 

1 

f = qT ltL (u L ) I] (Ti,l-iM) v .V t (26) 

j=L-l 

By using these decompositions, it is sufficient to keep in computer memory only 
four matrices: T 1L , T 1L _ 1 , M and V, with at most two nonzero elements per row, 
which strongly reduces amount of data to be stored. 

Since we want to get information about LR correlations, we should also take care 
of the connectivity between sites in the last column and the sites in the first column. 
In order to do so, one can imagine the top site in the last column of Z N to be some 
kind of "ghost" site representing the connection between the first site of last column 
of %n+l and the sites of the first column. Thus we need a column of height L to 
describe all connectivities in the model with the range of interaction equal to L-l. 
In order to keep the "ghost" site fixed, the matrix of the translation operator V will 
be constructed by acting with operator V on non-" ghost" sites only. Consequently, 
if we perform the following set of transformations 

U L ->• U L -1, -> U L -2, ■ ■ ■ , u 2 -> U!, — >• 0, (27) 

we will obtain ~T for the model with range L-l, with LR correlations included. 
With this redefinition, ~T has the following structure 

f=r- x 

V T + 

where ~T + and ~T_ contain the largest pL 1:L and the second largest pL 2 ,L eigenvalue 
of the whole matrix ~T, respectively. The correlation length is given by the relation 
(16) again. 
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The order of T for a model with range L is given by Nl, and the order of T + 
is given by N i _ 1 , where 

L+l 

N L =Y: V- (28) 
1=1 

= tf' 1 + <ffc* (29) 

for L>2, and (fl=3, (p r { = </?™ +1 = I. Some values of N L are presented in table 1. 

The present calculations were performed for the maximum range L=8 with 
q=l/16, 1/2, 1, 8, 16, 32 and 64. The largest eigenvalues of the matrices ~T±, 
are extracted by using the method explained earlier (section 2.2) with a similar 
accuracy and number of iterations. 

3 Results 

3.1 Convergence and extrapolation procedures 

The critical temperature and exponent u, given by the equations (7) and (8), still de- 
pend on the range M and it is important to apply the right extrapolation procedure. 
A simple analysis done within FSS (Privman and Fisher, 1983) can be reproduced 
in the case of the FRS (Glumac and Uzelac, 1989) and shows that the convergence 
is dominantly governed by the leading irrelevant field and the corresponding critical 
exponent y% < 0. On the basis of such analysis, we expect power-law convergence 
of T c and v in the large-M limit 

T CiM = T c + a- M^- 1 - (30) 

v- 1 = jy- 1 + b{T cM -T c )-M» + c • M yz = v- 1 + b' ■ M yz (31) 

where a, b, b' and c are constants. Thus, we chose to fit the obtained results for 
K c M and v M to the form 

y M = y e + A ■ M~ Xy (32) 

in the least-squares approximation (LSA), where y e denotes the extrapolated quan- 
tity. The calculations are performed by taking the five largest values of M. 
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It should be noticed, however, that Privman and Fisher (1983) analysis does 
not apply to the MF region and, consequently, the extrapolation form (32) becomes 
arbitrary there. On the other hand, in our previous analysis for the Ising model 
(Glumac and Uzelac, 1989) it turns out that the extrapolation form 

Vm = B + A(¥^)°* (33) 

gives a better agreement with the exactly known results for v in the MF region then 
the function of the form (32). 

The error bars of the above extrapolations can be estimated only roughly, by 
looking for the remaining L dependence of the results. In tables (2) and (3) are 
presented results which last digit is modified under the change of maximal range 
from L-l to L. The more careful examination of the errors, shows that the size of 
the error estimate changes with q and a, being less then ten percents at the most of 
q and a region. Exception is the small q and a region where the error is estimated 
to be a few times larger. 

3.2 Critical temperature 

The inverse critical temperature K Cj m, defined by equation (7), is calculated with 
an accuracy higher than 1CT 10 for different values of Potts states and with a as an 
parameter. 

When increasing the exponent of LR interaction, there will be a value of a 
beyond which the SR critical behaviour (with T c = 0) takes place. For the Ising 
model (q=2) it was analytically shown (Dyson, 1969) that this a is equal to 1. As 
argued more generally (Fisher et al. 1972, Sak 1973) this exchange of regimes should 
occur when the correlation function exponent i] of the LR system becomes smaller 
or equal to the SR one, i] SR . Thus in present case, we expect that a c = 1 for all 
values of q. This expectation is confirmed by our numerical results. In our earlier 
works (Glumac and Uzelac, 1989, 1991), we have been able to detect the appearance 
of SR forces governed critical behaviour by the change of K c> m from descending to 
ascending sequence and by the sudden increase of convergence exponent. In the 
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present calculations, both of these phenomena were observed for any q considered 
at a ~ 1. The extrapolated values for the critical temperature are presented in 
figures (la) and (lb). 

We have used the LSA extrapolation method with the form (32) to calculate 
K Cj€ and the leading convergence exponent Xk as a function of q and a which are 
presented in tables (2a) and (2b) respectively. The only point where this could not 
be applied is a = 1, where the simple convergence expression does not apply any 
more. In order to avoid non-monotonic behaviour, we have used only the data for 
the three largest values of M fitting a linear relation to them (x K = 1). The error 
should not be important, since for a = 1 the variation of data with M is rather 
weak. 

3.3 Critical exponent v 

The correlation length critical exponents vm has been calculated by using the equa- 
tion (8). The accuracy of results is reduced (with respect to that of the critical 
temperature) to the order of 1CT 6 due to the numerical differentiation. 

Following the form (31) for the scaling correction in non-MF region, one can 
expect that the correction terms to v^} would be smaller while calculating the v M 
on the temperature K ce instead of K c m- This is indeed the case for the q=2 
(s=l/2 Ising) model, where detailed comparison between vm(K c ^m) and vm(K c ^) 
has been given (Glumac and Uzelac, 1989). In the present model we observe better 
convergence of the data calculated from K Cfi in the LR a-region, for all q. 

According to the discussion at the end of section 3.1, the extrapolations in the 
low a region where MF behaviour is expected were performed at the temperature 
K c m using the extrapolation function of the form (33). Since the <Jmf border is still 
an open question for q > 2 (see later in text), we have used eq. (33) for a < .3, q < 1, 
and for a < .5, q = 2. 

The extrapolated results are presented in figure (2). In tables (3a) and (3b) are 
results for u^ 1 and the corresponding convergence exponents respectively. Similarly 
to preceding section, for a — 1 the LSA procedure was performed by imposing 
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x v = 1 in eq. (32). 

The special case q=2 which corresponds to the Ising model was extensively stud- 
ied in our earlier work (Glumac and Uzelac, 1989) with a maximum range equal 
10. The present work permits to reach the range of 20. Similar values of v e were 
obtained in both works, which suggests that an increase of range from 10 to 20 does 
not change the accuracy significantly. 

On table 4, our results are compared to the e-expansion results of Theumann 
and Gusmao (1985). A difference A of only few percent is obtained. Figure (2a) 
suggests that a = | as the MF border and a q-independent value of v^} F = a for all 
q< I- 

We have not found any analytical or numerical estimate for the values of o~mf 
or v for q > 3. In that situation we decided to calculate v at the temperature 
K ce since the scaling law ^ = const, is more accurate at that temperature than 
at K C ^ M . Also, we decided to make the extrapolation using function (32) which 
posesses a more transparent L-dependence. 

For large q (q > 16), only the continuous formalism could be applied, which 
does not permit to go beyond the range of L=8. Contrary to the small q case, for 
large q this range was not sufficient to find a good extrapolation. Namely, when 
applying the LSA there, with the assumed form (31), one obtains x v — > so that 
the correction terms become of the same order as the leading term and z/" 1 diverges. 
This behaviour of v suggests a change of regime of convergence and brings up the 
question of the order of transition. 

For all values of q the exponent v increases as a — > 1, which in the case q=2 
appeared (Glumac and Uzelac, 1989) as an indication of the essential singularity in 
this limit (Kosterlitz 1976). 

4 Conclusion and discussion 

By using a FRS method combined with transfer matrix calculations we have been 
able to study the long-range Potts model on the one-dimensional lattice in a wide 
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range of values for the number of states q and the interaction exponent a. We have 
shown that for this problem the transfer matrix can be decomposed into matrices 
of much simpler form. A different transfer matrix procedure can be done also in the 
case of continuous q, so that our analysis could be extended to some non-integer q 
cases of interest. 

The study is concentrated on the phase diagram and the correlation length crit- 
ical exponent v. 

The critical temperature shows monotonic decrease with q and a and has a finite 
value at a — 1. The change of behaviour of K c M at a ~ 1 indicates the exchange 
from LR to SR regime, established for any q. The critical exponent was calculated 
at the extrapolated critical temperature, and good agreement with values obtained 
by other authors was observed. Generally, a good accuracy is harder to obtain in 
the small-a region due to the very long-range of interaction and in the a — > 1 region 
where the range of interaction ceases to be a good scaling variable. 

There are two points that require some further discussion. 

First is a question of mean-field border. The present method, unlike the FSS, 
has the advantage that it can be applied within the mean-field region, but on the 
other hand, passes smoothly between the two regions (Glumac and Uzelac 1989, 
1991) , which does not permit to point out the mean-field border omf- Another 
way to detect this border could be numerical, by observing the change of behaviour 
of v as a function of a. But, our previous results on the Ising model where <jmf 
is known, show that close to it the numerical results for v are not sharp enough to 
locate this point with precision. In the present case, we can only confirm that our 
results are consistent with o MF as known from the literature for a few particular 
cases: q=2 [omf — \ Fisher et al 1972), q=l {omf = § Priest and Lubensky 1976), 
Q. — \ (%f = | Aharony 1978, Aharony and Pfeuty 1979), but the precise location 
of omf for arbitrary q is left unestimated. 

A second point is the possible appearance of a first order transition for some 
q > q c . Within the transfer matrix formalism this problem was considered for the 
2D SR case where q c exists and is exactly known (q c = 4, Baxter 1973). Igloi 
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and Solyom (1983) have shown that the first order transition is connected with the 
crossing of the largest and the third largest eigenvalue of the transfer matrix for 
finite length chain and large q. We have applied a similar calculation in our model 
considering two groups of parameters (q=100, a = .8, L=4 and q=300, a = .4, L=4). 
Both cases give a negative result; the first and third eigenvalues do not intersect as 
a function of temperature. By analogy with Igloi and Solyom, that result can be 
interpreted as the absence of a first-order transitions for any q in the ID LR Potts 
model. Since this is opposite to the indications of the behaviour of v for large q (at 
the end of section 3.3), the question of the order of transition for large q is still left 
open. 
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Figure captions: 
Figure la: 

The extrapolated critical temperature as a function of a. 
Figure lb: 

The extrapolated critical temperature as a function of a. Figures (a) and (b) 
have common q=l line which allows the comparison of scales on both figures, 
solid filled circle; 

Figure 2a: 

The exponent v~ l as a function of a. solid filled circle. 
Figure 2b: 

The exponent v~ x as a function of a. Figures (a) and (b) have common q=2 line 
which allows the comparison of scales on both figures. 
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Table captions: 



Table 1: 

N L is the order of transfer matrix for continuous-q model with maximum range 
equal L. 

Table 2a: 

The extrapolated values of inverse critical temperature as a function of q and a. 
For the error bars see in text. 

Table 2b: 

The critical temperature convergence exponent xk as a function of q and a. 
Table 3a: 

The extrapolated values of u^ 1 as a function of q and a. For the error bars see 
in text. 

Table 3b: 

The convergence exponent Xj; CIS 8b function of q and a. 
Table 4: 

The comparison between extrapolated values u^ 1 and Theumann and Gusmao 
values i>tg for a = .4. 
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